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Abstract — In this paper we report some extraordinary non- 
linear dynamics measured in the resonance curve of NbN 
superconducting stripline microwave resonators. Among the 
nonlinearities observed: aburpt bifurcations in the resonance 
response at relatively low input powers, asymmetric resonances, 
multiple jumps within the resonance band, resonance frequency 
drift, frequency hysteresis, hysteresis loops changing direction 
and critical coupling phenomenon. Weak links in the NbN grain 
structure are hypothesized as the source of the nonlinearities. 

Index Terms — Bifurcations, nonlinear effects, hysteresis, NbN, 
microwave resonators. 
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I. Introduction 

NONLINEAR effects in superconductors in the microwave 
regime have been the subject of a large number of 
intensive studies in recent years. Most of the attention is 
focused on studying one or more of the following issues: in- 
vestigating the origins of nonlinear effects in superconductors 
[1],[2], introducing theoretical models that explain nonlinear 
behavior [3], [4], identifying the dominant factors that mani- 
fest these effects [5], [6], find ways to control and minimize 
nonlinear effects [7], [8] such as, harmonic generation and 
intermodulation distortions, which degrade the performance 
of promising superconducting microwave applications mainly 
in the telecommunication area [9]. 

Among the nonlinear effects reported in the literature as- 
sociated with resonance curves, one can find the commonly 
known Duffing nonlinearity which is characterized by skewed 
resonance curves above certain power level, appearance of 
infinite slope in the resonance lineshape, pronounced shift of 
the resonance frequency and hysteretic behavior [10],[1 1],[12]. 
To account for this effect, associated with the rise of kinetic 
inductance of superconductors, both thermal [13], and weak 
link [10] explanations have been successfully applied. Other 
nonlinear effects were reported by Portis et al. [14], where 
they observed notches that develop on both sides of the 
frequency response of their HTS microstrip patch antenna, 
accompanied with hysteresis and frequency shift, as they 
have driven their antenna into the nonlinear regime. Similar 
results were reported also by Hedges et al. [15], in their 
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YBCO stripline resonator, and by [16] in their YBCO thin 
film dielectric cavity. All three studies attributed the observed 
nonlinear behavior to abrupt changes in the resistive loss of 
weak links, thermal quenching and weak link switching to 
normal state. 

In this study, being interested in the behavior of nonlinear 
resonances, we have fabricated different NbN superconduct- 
ing microwave resonators exhibiting some unusual nonlinear 
effects, which to the best of our knowledge, have not been 
reported before in the literature. We study the dependence of 
these resonators on the injected input power level, and examine 
the resonance curve behavior under different scan directions 
showing interesting features. To account for our results, we 
consider briefly some possible physical mechanisms that may 
be responsible for the observed effects. 

II. Resonators Design 

A. Resonator Geometries 

The resonators were designed in the standard stripline 
geometry, which consists of five layers as shown in the cross 
section illustration depicted in Fig. \l\ 

The superconducting resonator was dc-magnetron sputtered 
on one of the sapphire substrates, whereas the superconducting 
ground planes were sputtered on the inner covers of a gold 
plated Oxygen Free High Conductivity (OFHC) Cupper Fara- 
day package that was employed to house the resonators. The 
dimensions of the sapphire substrates were 34mm X 30mm 
X 1mm. The resonator geometries implemented, which we 
will refer to them, for simplicity, by the names Bl, B2, B3, 
are presented in the insets of Figs. [3] |4] |5] respectively. The 
width of the feedlines and the thin part of the resonators 
was set to 0.4mm to obtain characteristic impedance of 5057. 
The gap between the feedline and the resonators was set to 
0.4mm in the Bl, B3 cases, and to 0.5mm in B2 case. The 
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TABLE I 
Sputtering Parameters 



Process parameter 


Bl 


B2 


B3 


Partial flow ratios (Ar,N2) 


(87.5%,12.5%) 


(75%,25%) 


(70%,30%) 


Base temperature 


11°C 


11 °C 


13°C 


Total pressure 


6.9 ■ 10~ 3 torr 


8.1 • 10" 3 torr 


5.7 • 10" 3 torr 


Discharge current 


0.36A 


0.55A 


0.36A 


Discharge voltage 


351V 


348V 


348V 


Discharge power 


121W 


185W 


133W 


Deposition rate 


6^ 

sec 


7.8^ 

sec 


3.8^ 

sec 


Thickness (t) 


2200A 


3000A 


2000A 


Base pressure 


3.1 ■ 10" 8 torr 


7.3 ■ 10" 8 torr 


8 ■ 10 -8 torr 


Target-substrate distance 


80mm 


90mm 


90mm 



frequency modes of Bl, B2, B3 resonators were theoretically 
calculated using a simple transmission line model, presented 
in appendix A, and were also experimentally measured using 
vector network analyzer (NA). The theoretical calculation was 
generally found to be in good agreement with the measurement 
results, as discussed in Appendix A. 

III. Fabrication Process 

The sputtering of the NbN films was done using a dc- 
magnetron sputtering system. All of the resonators reported 
here were deposited near room temperature [17], [18], [19], 
where no external heating was applied. The system was usually 
pumped down prior to sputtering to 3 — 8 • 10 _8 torr base 
pressure (achieved overnight). The sputtering was done in 
Ar/N2 atmosphere under current stabilization condition [20]. 
The relative flow ratio of the two gases into the chamber and 
the total pressure of the mixture were controlled by mass flow 
meters. The sputtering usually started with a two minute pre- 
sputtering in the selected ambient before removal of the shutter 
and deposition on the substrate. The sputtering parameters of 
the three resonators are summarized in table 1. Following the 
NbN deposition, the resonator features were patterned using 
standard photolithography process, whereas the NbN etching 
was done using Ar ion-milling. 

The fabricated resonators were characterized by relatively 
low T c for NbN and relatively high normal resistivity p, that is 
in good agreement with Ref. [18]. T c measured for Bl, B2 and 
B3 was 10.7, 6.8, 8.9 [K] respectively, whereas p measured 
for B2 and B3 was 348 and 500 [/iSlcm] respectively. 

To obtain resonators with reproducible physical properties 
we have used the sputtering method discussed in [20], where it 
was claimed that reproducible parameters of films are assured, 
by keeping the difference between the discharge voltage in 
a gas mixture, and in pure argon, constant, for the same 
discharge current. In Fig. |2] we show one of the characteri- 
zation measurements applied to our dc-magnetron sputtering 
system, exhibiting a knee-shape graph of discharge current as 
a function of discharge voltage. The knee-shape graph was 
obtained for different Ar/N2 mixtures at room temperature as 
shown in the figure. The discharge voltage difference measured 
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Fig. 2. Discharge current vs. discharge voltage of the sputtering system dis- 
playing current- voltage knee for different percentages of Argon and Nitrogen, 
at ambient temperature of 20° C. 



in the presence of N2 gas relative to the value measured in pure 
argon at the same discharge current, is also pointed out in the 
figure, corresponding to different currents and N2 percentages. 

IV. Measurement Results 

All measurements presented in this paper have been con- 
ducted at liquid helium temperature. 

A. S\\ Measurements 

The resonance response of the resonators was measured 
using the reflection parameter Snof a vector NA. The res- 
onance response obtained for the third mode of Bl resonator 
~8.26GHz at low input powers, between —23 dBm and —18 
dBm in steps of 0.05 dBm, is shown in Fig. [3] A small offset 
was applied between the sequential graphs, corresponding 
to different input powers, for clarity, and to emphasize the 
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Fig. 3. Sn measurement of Bl resonance at ~ 8.26GHz with 10MHz 
span, at low input powers, exhibiting extraordinary nonlinear effects. The 
resonance curves corresponding to different input powers were shifted by a 
constant offset for clarity. The inset shows the resonator geometry. 



Fig. 4. Sn measurement of B2 resonance at ~ 4.385GHz with 30MHz 
span, at low input powers, exhibiting strong nonlinear effects. The resonance 
curves corresponding to different input powers were shifted by a constant 
offset for clarity. The inset shows the resonator geometry. 



nonlinear evolution of the resonance response as the input 
power is increased. The interesting characteristics of this 
nonlinear evolution could be summarized as follows: 

1) In the power range between —23.5 dBm and —23.25 
dBm, the resonance is symmetrical and broad. 

2) At input power level of —23.25 dBm, a sudden jump of 
about — 15 dB occurs in the resonance curve at the minima 
where the slope of the resonance response is small. 

3) As the input power is increased in steps of 0.05 dBm 
the resonance becomes asymmetrical, and the left jump shifts 
towards the lower frequencies gradually. 

4) As we continue to increase the input power, the jumps 
decrease their height but the resonance curve following the 
jumps becomes more symmetrical and deeper, and at certain 
input power level we even witness a critical coupling phe- 
nomenon where Su (oj) at resonance is almost zero, no power 
reflection is present. 

5) The resonance becomes symmetrical again and broader 
and the bifurcations disappear. 

6) All previously listed effects occur within a frequency 
span of 10MHz, power range of about 5 dBm, and power step 
of 0.05 dBm. 

Moreover, in order to estimate how narrow this resonance 
could be in the vicinity of critical coupling, a separate Su 
measurement has been applied near critical coupling power, 
using 1601 measurement points and frequency span of 2 MHz. 
The bandwidth of the resonance curves measured, according 
to the +3 dB method, was about 0.25 • 10~ 5 GHz, whereas the 
ratio //A/~ 3.3- 10 6 . 

Similar behavior to that exhibited by the nonlinear third 
mode of resonator Bl can be clearly seen in Fig. 0] and Fig. 
[5J which show the nonlinear dynamic evolution of the second 
mode of resonator B2 and the first of B3 respectively. The 
main differences between the figures are: 

1) The power levels at which these nonlinear effects take 
place. Whereas in Bl case they happen between —23 dBm 




1.4 1.45 1.5 1.55 1.6 1.65 1.7 1.75 1.8 
Frequency [GHz] 



Fig. 5. Nonlinear response of B3 resonance at ~ 1.6 GHz, corresponding to 
input power levels increasing in 0.01 dBm. At input power of 1.49 dBm we 
observe two obvious bifurcations in the resonance band, and another small 
one at the right side, marked with circles. The different resonance curves 
were shifted by a constant offset for clarity. The inset shows the resonator 
geometry. 



and —18 dBm, in B2 case they happen between —9.5 dBm 
and —1.5 dBm and in B3 case they happen around 1 dBm. 

2) In Fig. [5] corresponding to B3 resonator we witness two 
apparent bifurcations within the resonance band as indicated 
by circles on the figure, a feature that we did not encounter 
in Figs. |3] and 0] 

B. Verifications 

In order to verify that the bifurcation feature, previously 
measured using NA Su parameter, is not a measurement 
artifact, we applied a different measurement configuration, 
shown in Fig. |S] where we scanned the frequency axis with 
CW mode of an Agilent synthesizer and measured the reflected 
power from the resonator by a power diode and voltage meter. 
The load that appears in Fig. [6] following the diode is an 
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Fig. 6. A setup that was used to verify the occurrence of the bifuractions at 
the nonlinear resonance (Bl third mode), measured previously by NA. 
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Fig. 7. Nonlinear response of Bl third mode measured using the configuration 
shown in Fig. |6| The red line represents a forward scan whereas the blue line 
represents a backward scan. Two abrupt jumps appear at both sides of the 
resonance curve and small hysteresis loops are present at the vicinity of the 
jumps. The resonance curves corresponding to different input powers were 
shifted by a constant offset for clarity. 



Agilent load used in order to extend the linear regime of the 
power diode. The results of this measurement configuration 
are shown in Fig.0 The frequency scan around the resonance 
was done using 201 points in each direction (forward and 
backward). A small hysteresis loop can be seen around the 
two bifurcations. 



C. Abrupt Bifurcations 

In attempt to find out whether the resonance curve of these 
nonlinear resonances changes its form along two or more 
frequency points, further measurements where carried out 
using NA, where we scanned the frequency axis in the vicinity 
of the jump with high frequency resolution. The measurement 
results corresponding to frequency step of ^600 Hz and ~2.5 
kHz are presented in Fig. [8] plot (a) and (b) respectively, 
indicating abrupt transition between two bistable states. 



Fig. 8. (a) A forward CW mode scan using NA, in a 12MHz span around 
the left jump of the nonlinear resonance of B2. The scan includes 20,000 
frequencies, which is equivalent to a frequency step of ~600Hz between the 
data points. In spite of this small frequency step, the jump still occurs between 
just two sequential frequencies, (b) A backward CW mode scan using NA, in 
a 10MHz span around the left jump of the nonlinear resonance of B2. The 
scan includes 4000 frequencies, which is equivalent to a frequency step of 
~2.5kHz between the data points. Also for this case the jump occurs between 
just two sequential frequencies. 



D. Frequency Hysteresis 

Applying forward and backward frequency sweeps to these 
resonators, reveals a very interesting hysteretic behavior. In 
Fig. [9] we show a representative frequency scan of B2 second 
mode, applied in both directions, featuring the following 
nonlinear dynamic behavior: 

1) At low input powers —8.05 dBm and —8.04 dBm, the 
resonance is symmetrical and there is no hysteresis. 

2) As the power is increased by 0.01 dBm to —8.03 dBm, 
two bifurcations occur at both sides of the resonance response 
and hysteresis loops form at the bistable regions. 

3) As we continue to increase the input power gradu- 
ally, the hysteresis loop, associated with the right bifurca- 
tion, changes direction. At first it circulates counterclockwise 
between -8.03 dBm and -7.99 dBm, at -7.98 dBm the 
two opposed bifurcations, at the right side, meet and no 
hysteresis is detected, as we increase the power further the 
right hysteresis loop appears again, circulating, this time, in 
the opposite direction, clockwise. 

Furthermore, it is worth mentioning that the hysteresis loops 
changing direction are not unique to this resonator, or to the 
bifurcation occurring on the right side of the resonance. It 
appears also in the modes of Bl, and it occurs at the left side 
bifurcation as well, but at different power level. 

E. Multiple Bifurcations 

Frequency sweep applied to B3 first resonance in both 
directions, exhibits yet another feature, in addition to the two 
bifurcations at the sides of the resonance curve, which we 
have seen earlier, there are another two smaller bifurcations 
accompanied with hysteresis within the resonance lineshape, 
adding up to 4 bifurcations in each scan direction, as can be 
seen in Fig. This feature may have a special significance 
in explaining the physical origin of these nonlinearities. 
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Fig. 9. Forward and backward scan measurement, performed using NA, 
measuring B2 second mode nonlinear resonance. The red line represents a 
forward scan whereas the blue line represents a backward scan. The graphs 
exhibit clear hysteresis loops forming at the vicinity of the bifurcations, 
and hysteresis loop changing direction as the input power is increased. The 
resonance curves corresponding to different input powers were shifted by a 
constant offset for clarity. 
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Fig. 10. Sii parameter measurement of the first resonance of B3 at input 
powers 1.52 dBm through 1.58 dBm in steps of 0.02 dBm using forward 
and backward CW mode scan of NA employing 2000 measurement points 
in each direction. The graph shows clearly 4 jumps within the band of the 
resonance in each direction. The resonance curves corresponding to different 
input powers were shifted by a constant offset for clarity. 



V. Discussion 

This unusual dynamic behavior of our NbN films, demon- 
strated earlier, highly suggest built-in Josephson junctions, 
forming at the boundaries of the granular NbN columnar 
structure [21], [22], as the underlying physical mechanism re- 
sponsible for the observed nonlinearities [1]. Another physical 
mechanism that may be considered as a strong candidate for 
explaining the effects, is the local heating mechanism which 
was hypothesized as the source of notches and switching 
effects observed in HTS films [14], [15], [16]. Nevertheless 
recent measurements done on SQUID ring containing a single 
Josephson junction inductively coupled to a radio frequency 
resonant circuit [23], [24], [25] exhibiting opposed bifurca- 
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tions, hammerhead resonance lineshape and similar effects, 
suggest similar physical mechanism and further support the 
former hypothesis. 

VI. Conclusion 

In the course of this experimental work we have fabricated 
several stripline NbN resonators dc-magnetron sputtered on 
sapphire substrates at room temperature implementing dif- 
ferent geometries. The resonators have exhibited similar and 
unusual nonlinear effects in their resonance response curves. 
The onset of the nonlinear effects in these NbN resonators 
varied between the different resonators, but usually occurred 
at relatively low powers, typically 2-3 orders of magnitude 
lower than Nb for example. Among the nonlinear effects 
observed: abrupt and multiple bifurcations in the resonance 
curve, power dependent resonance frequency shift, hysteresis 
loops in the vicinity of the bifurcations, hysteresis loops 
changing direction, and critical coupling phenomenon. Weak 
links forming in the NbN films are hypothesized as the source 
of the nonlinearities. Further study of these effects under other 
modes of operation and measurement conditions would be 
carried in the future, in order to substantiate our understanding 
of these extraordinary effects. 

Appendix I 

Resonance frequency calculation of Bl, B2, B3 

RESONATORS 

The calculation process of the resonance frequencies of Bl 
and B2 makes use of opposite traveling voltage-current waves 
method [26], [27]. For this purpose we model Bl and B2 
resonators as a straight transmission line extending in the z- 
direction with two characteristic impedance regions Z\ and 
Z-i as shown in Fig. \H\ 

The equivalent voltage along the resonator transmission line 
would be given, in general, by a standing waves expression in 
the form: 

!A + cos f3 (z — a) + B + sin j3{z — a) ze(a, b) 
A cos [3z + B sin [3z ze(—a,a) 
A- cos f3(z + a) + £>_ sin (3 (z + a) ze(—b, —a) 

_ (1) 
where (3 (= 27r/y / e7/c) is the propagation constant along 

the transmission line, and A + , B + , A, B, A_, B_ are con- 
stants that can be determined using boundary conditions and 
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applied power amplitude. However due to the symmetry of 
the problem z < — ► —z, we expect the solutions to have 
defined parity, where V (z) — V (— z) for symmetric solution 
and V (z) — —V (— z) for antisymmetric solution. Thus by 
taking advantage of this property and demanding that V (z) 
be continuous at z = a and z = —a, one gets: 
cos Pa cos P (z — a) 
+B S sin/3 (z — a) 

V S yrn (z) = { COS (3z 

cos Pa cos P (z + a) 
—B s sin/3 [z + a) 
sin Pa cos P (z — a) 
+B a sin/3 (z — a) 
Vanti (z) = { sin j3z 

— sin Pa cos P (z + a) 
+B a sin/3 (z + a) 



ze(a, b) 
ze(—a, a) 
ze(—b, —a) 

ze(a, b) 
ze(—a, a) 
ze(—b, —a) 



where V sym (z) stands for the symmetric solution whereas 
Vanti (z) for the antisymmetric solution. To calculate the value 
of the new constants B s ,B a , we require that the equivalent 
current / (z) along the transmission line, which is given 
by I (z) — (ij pZi)dV I dz where Z% is the characteristic 
impedance of the line, be continuous at z = a and z = —a. 
Following this requirement one gets B s = —i] sin (/3a) and 
B a = V cos (Pa), where r] = Z2/Z1. The symmetric and 
antisymmetric solutions of V (z) and I (z) are given by: 
cos Pa cos P(z — a) 
-r] sin Pa sin /3 (z — a) 

Vsyrn (z) = { COS Pz 

cos Pa cos P (z + a) 
+rj sin Pa sin /3 (z + a) 
sin Pa cos P(z — a) 
+ry cos Pa sin P (z — a) 
Vanti {z) = { sin (pz) 

— sin (Pa) cos P (z + a) 
+i] cos (Pa) sin P(z + a) 

-i cos /3a sin f3(z — a) 



ze(a, b) 
ze(—a, a) 
ze(—b, —a) 

ze(a, b) 
ze(—a, a) 
ze(—b, —a) 



' s y rn 



(z) 



( z 2 

in sin (3a cos (3(z — a) \ 

Bigg) 



(*)=< 



Zi 

/ — i cos /3a sin (3(z+a) 

v zi 

. in sin (3a cos /3(z+a) \ 

H Z~ 2 > 

/ —2 sin (3a sin (3(z — a) 

, ir; cos /3a cos 0(z— a) « 

Z2 
i cos (3z 

Zx 

/ i sin /3a sin /3(z+a) 

. in cos (3a cos /3(2;+a) ^ 
H Z 2 - 



ze(— a, a) 
ze(— 6, —a) 

ze(a, b) 
ze(—a, a) 
ze(—b, —a) 



A. Bl Resonator 

Since the resonator ends are shorted we demand V (b) = 
V(-b). 

The symmetric case: 

In this case we either have maximum or minimum at V (b) 
thus we get / (b) = 0, yielding the following condition on the 
resonance frequencies: 

cos (Pa) sin (/3 (b — a) ) + r\ sin (/3a) cos (P (b — a) ) = (2) 



The antisymmetric case: 

From the antisymmetric V (— b) — — V (— b) and the conti- 
nuity V (b) = V (-6) conditions, we get V (b) = V (-6) = 
0, which yields: 



sin (Pa) cos (P(b— a)) + 77 cos (/3a) sin (/3 (b — a)) = (3) 

Substituting the following numerical values ?/ = Z^jZi = 
0.5, Zi = a = 13mm, I2 = b — a = 6.5mm into 
the above resonance frequency conditions and solving for 
frequencies below 10GHz, yields the following solutions 
(2.5035GHz, 5.697GHz, 8.1647GHz) to the symmetric case, 
and (2.9804GHz, 5.1786GHz, 8.1647GHz) to the antisym- 
metric case, with doubly degenerate mode at 8.1647GHz. By 
comparing these calculated resonances to the directly mea- 
sured resonances of B 1 resonator, obtained using a broadband 
Sn measurement (2. 5812GHz, 5.6304GHz, 8.4188GHz), we 
find that the excited resonances correspond to the symmetrical 
case only. The antisymmetric modes do not get excited because 
they have a voltage node at the feeding line position. 

B. B2 Resonator 

Since the resonator ends are open-circuited we demand 
I(b) = I(-b)=0. 
The symmetric case: 

We require that the current associated with the symmetric 
voltage, vanishes: 

cos (Pa) sin \fl (b - a)} + r\ sin (/3a) cos [/3 (b - a)} = (4) 
The antisymmetric case: 

We require that the current associated with the antisymmet- 
ric voltage, vanishes: 

— sin (Pa) sin (ft (b — a)) + r] cos (Pa) cos P (b — a) — (5) 

Substituting the following numerical values rj = Z2/Z1 = 
49.9/10.4 = 4. 79, h = a = 11.97mm, l 2 = b - a = 6.43mm 
into the above resonance frequency conditions and solving 
for frequencies below 10GHz, yields the following solutions 
(2.6486GHz, 6.0288GHz, 8. 5588GHz) to the symmetric case, 
and (1.1763GHz, 4.3698GHz, 7.4597GHz, 9. 7778GHz) to 
the antisymmetric case. By comparing these calculated reso- 
nances to the directly measured resonances of B2 resonator, 
obtained using a broadband Sn measurement (2. 5152GHz, 
4.425GHz, 6.3806GHz, 8.176GHz), we find a good agree- 
ment between the two results. The missing resonances do not 
get excited apparently because of the coupling location of the 
feedline relative to the resonator. 

C. B3 resonator 

B3 resonator, in contrast, showed some larger discrepancy 
between the measured value for the first mode ~ 1.6 GHz 
(seen in Fig. and the theoretical value f\ = 2.4462 GHz 
calculated according to the approximated equation: 

nc 



In 



2/ £ y 



(6) 
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where n is the mode number, c is the light velocity, I is the 
open-circuited line length (~ 20mm), and e r is the relative 
dielectric coefficient of the sapphire (~ 9.4). 
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